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The Oubangui is a major tributary of the Congo River. We describe the biogeochemistry of contrasting 
tributaries within its central catchment, with watershed vegetation ranging from wooded savannahs to 
humid rainforest. Compared to a 2-year monitoring record on the mainstem Oubangui, these tributaries 
show a wide range of biogeochemical signatures, from highly diluted blackwaters (low turbidity, pH, 
conductivity, and total alkalinity) in rainforests to those more typical for savannah systems. Spectral 
analyses of chromophoric dissolved organic matter showed wide temporal variations in the Oubangui 
compared to spatio-temporal variations in the tributaries, and confirm that different pools of dissolved 
organic carbon are mobilized during different hydrological stages. 5 13 C of dissolved inorganic carbon 
ranged between — 28.1%o and — 5.8%o, and was strongly correlated to both partial pressure of C0 2 and to the 
estimated contribution of carbonate weathering to total alkalinity, suggesting an important control of the 
weathering regime on C0 2 fluxes. All tributaries were oversaturated in dissolved greenhouse gases (CH 4 , 
N 2 0, C0 2 ), with highest levels in rivers draining rainforest. The high diversity observed underscores the 
importance of sampling that covers the variability in subcatchment characteristics, to improve our 
understanding of biogeochemical cycling in the Congo Basin. 



The recognition that carbon (C) processing within inland waters could be a substantial component in C 
budgets at the catchment, regional, or global scale (e.g. 13 ), has led to an increased momentum in studies on 
riverine biogeochemistry. Further constraining the role of river systems will only be possible by generating 
additional comprehensive datasets, in particular on systems or regions which are currently underrepresented. Of 
key importance are tropical and subtropical regions, which have been suggested to be of particular importance in 
terms of riverine transport of sediments and carbon (e.g. 4 5 ), and have been suggested to show higher areal C0 2 
outgassing rates than their temperate or boreal counterparts 3,6 . Based on a long tradition of interdisciplinary 
studies, the Amazon often serves as a model system for the biogeochemical functioning of tropical river basins 
with African rivers being largely neglected. The Congo River basin is second to the Amazon in terms of discharge 
(1310 versus 6640 km 3 yr" 1 ) and catchment size (3.8 10 6 versus 5.8 10 6 km 2 ), but is much less well characterized 
from a hydrological, physico-chemical, and biogeochemical point of view. The main biogeochemical studies on 
the Congo basin have focussed on dissolved and particulate erosion rates both on the main Congo River and in 
some of its major tributaries such as the Oubangui and Sangha (e.g., 711 ). The synthesis of Laraque et al. 11 , however, 
shows that sampling efforts were concentrated on a limited number of sites, and these monitoring programmes 
were limited to a restricted number of basic parameters. For many stations, extrapolations needed to be made 
based on less than 5 samples (e.g. the upper Congo River and Kasai River 12 ,). Mario tti et al. 13 provided the first 
stable isotope data on particulate organic carbon (POC) from the lower Congo River and some of its tributaries. 
More recently, some explorative work on the organic and inorganic C biogeochemistry of the lower Congo River 
has been presented 14 " 46 , and we reported results from a 1-year sampling programme on the Oubangui at Bangui 
(Central African Republic, CAR), focussing on the dynamics and annual export of different C pools and green- 
house gas (GHG) emissions 17 . Rivers in the Congo basin network, however, can be expected to show highly 
variable physico-chemical and biogeochemical characteristics, considering the range in geology, climatic condi- 
tions and land use/vegetation across the basin. Laraque et al. 18 analysed the physico-chemical properties and major 
ion composition in different tributaries of the lower Congo, and found strong contrasts between those draining the 
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"Cuvette Congolaise" (mainly rainforest- dominated lowlands) and 
the "Teke plateaux", where savannahs form the main biome type. 
Similarly, Mann et al. 19 recently demonstrate strong gradients in 
biogeochemical characteristics across rivers and streams in the 
Sangha and lower Oubangui catchment. The existence of a large 
gradient in biogeochemical characteristics is relatively well documen- 
ted for some of the large South American basins (e.g. 20 ), but compre- 
hensive geochemical data from the Congo basin are still mostly 
limited to the mainstem in the Malebo Pool at Brazzaville and some 
of the main tributaries. During three sampling campaigns between 
2010 and 2012, we sampled a number of tributaries of the Oubangui 
River, from savannah- dominated systems (Mbali River, Mpoko 
River) to those draining mainly humid rainforest ecosystems 
(Lobaye, Mbaere, Bodingue, and some of their minor tributaries) in 
the Ngotto Forest (Figure 1). Here, we compare biogeochemical data 
from these tributaries with our own results from 2 years of high- 
frequency sampling on the mainstem Oubangui River at Bangui, and 
with available literature data on the wider Congo Basin. 

Results 

Material fluxes in the Oubangui River. A detailed discussion of 
results for the first year of monitoring was presented in Bouillon 
et al. 17 . While no direct DOC measurements are available for the 
2 nd year, we estimated DOC fluxes based on the relationship between 
a 35 o and DOC observed in the tributaries (see Supplement Figure SI; 
unpublished data from other sites within the larger Congo Basin also 
follow this trend). The 2 nd year of measurements was characterized 
by a markedly lower annual discharge (—30% lower), and for most 
elements, the annual transport fluxes were accordingly lower than 



those reported previously by 22% (for DIC) to 55% (for DOC, 
summarized in Table 1). As an illustration, the seasonal patterns 
over both years of sampling are presented in Figures 2-3 for TSM, 
TA, 5 13 C DIC , pC0 2 , CH 4 and N 2 0. The full data can be found in 17 
and in Supplementary Table 2. During the 2 nd year, however, we also 
analysed spectral parameters of cDOM. The a 350 values exhibited a 
strong increase with increasing discharge, with values ranging from 
about 5 m" 1 during base flow to > 20 m" 1 during high flow 
(Figure 4). The a 2 5o:a 365 ratio values were relatively high during 
base flow (4.6-5.5) and dropped markedly during the ascending 
limb of the hydrograph to values closed to 4.2 during high flow 
period (Figure 4). Opposite patterns were observed for the spectral 
slopes. S 275 _295 values decreased with increasing discharge (from 
>0.016 nm -1 during base flow to 0.012 nm" 1 during high flow) 
whereas S 350 _4oo values increased with increasing discharge (from 
0.013 nm" 1 during base flow to 0.016 nm" 1 during high flow). 
Finally, spectral slope ratios (S R ) exhibited a strong seasonal 
variation with markedly higher values during base flow period 
(>1.2) compared to high flow period (about 0.8). 

Despite the lower transport fluxes, estimated C0 2 fluxes across the 
water-air interface were similar or slightly higher during 2011-2012 
than during 2010-201 1 (Table 1). Overall, the seasonality in concen- 
trations and isotope ratios followed similar patterns during both 
years (Figures 2-3), hence we will not discuss these in detail but 
use the data mainly to examine overall patterns across the mainstem 
and tributaries. 

Tributary characteristics. Tributaries showed a much wider range 
of values for most of the measured variables than recorded during the 
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Figure 1 | Map showing the location of the Oubangui within the Congo Basin, and detailed map showing the tributary sampling sites (note that some 
sites represent small first order tributaries). Map produced using ArcGIS software. 
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Table 1 | Summary of fluxes and element ratios, and yields (i.e. fluxes expressed per unit are of the catchment) measured in the Oubangui 
River at Bangui during two subsequent hydrological years, 20th March 201 0 to 1 9th March 201 1 (data from Bouillon et al. 201 2) and 20th 
March 201 1 to 1 8th March 201 2 (this study). Fluxes are expressed in Tg y" 1 , yields are expressed in ton km" 2 y" 1 

Constituent Flux, 201 0-201 1 Flux, 201 1 -201 2 Yield, 201 0-201 1 Yield, 201 1 -201 2 



Discharge (m 3 s ] ) 


3361 


2343 






TSM 


2.45 


1.63 


5.01 


3.33 


POC 


0.1446 


0.0889 


0.296 


0.182 


PN 


0.0142 


0.0100 


0.029 


0.020 


DOC 


0.725 


0.323 a 


1.483 


0.661 


DIC 


0.479 


0.374 


0.980 


0.765 


%POC 


5.9 


5.5 






POC/PN (mass based) 


10.2 


8.9 






Annual water-air CO2 flux (mmol CO2 m" 2 d" 1 ) 1 


24.5 b -46.7 c 


28.9 b -53.3 c 







a: estimated based on relationship between 0350 and DOC concentrations, see text and Supplementary Figure SI . b ; c: calculated using two different parameterisations, see 17 for details. 



2 years of monitoring on the mainstem Oubangui River. 
Biogeochemical signatures in the tributaries appear to be clustered 
according to the catchment characteristics. Specific conductivity 
generally showed very low values (average 14.6 |iS cm" 1 , with a 
minimum of 7.6 jaS cm" 1 ) in many sites within the densely 
vegetated Ngotto Forest (Lobaye, Mbaere, Bodingue), and was 
higher in the more northern, savannah catchments (Mbali, Mpoko 
: 22-170 |iS cm" 1 ) as well as in all left bank tributaries of the Lobaye 
river (average 89.4 |iS cm" 1 ) which are situated at the rainforest/ 
savannah transition. TA was generally well correlated with specific 
conductivity (Figure 5), and ranged between 0.009 and 1.791 mmol 
L" 1 . The highly dilute rivers in rainforest- dominated catchments 
were also characterised by very low TSM concentrations (0.7- 
16.0 mg L" 1 ) while higher TSM loads were encountered in some of 
the tributaries of the Lobaye (9.4-44.4 mg L" 1 ) and in Mpoko River 
(78.0 mg L" 1 ). Rivers in the Ngotto Forest were also characterized by 
relatively lower pH values, but despite the visual appearance as 
'blackwater' rivers (in particular the Mbaere and Bodingue), do not 
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Figure 2 | Seasonal variations of daily discharge (dotted lines) and (A) 
total suspended matter concentrations (TSM), and (B) concentrations of 
CH 4 and N 2 0 in the Oubangui River at Bangui between March 2010 and 
March 2012. 



show substantially higher DOC concentrations (between 1 and 6 mg 
L" 1 for the majority of samples, Figure 6). The a 350 values of the 
Oubangui tributaries were higher during the wet period (36.2 ± 
14.0 m" 1 ) than during the dry period (22.9 ± 6.8 m" 1 ) (see 
Supplementary Tables) and are well correlated with DOC 
concentrations across the seasons (Figure SI). The variations of 
a 250 :a 365 ratios were relatively limited within and across seasons, 
with a pattern of lower values during wet period (3.84 ± 0.16) 
compared to the dry period (4.02 ± 0.17). Spectral slopes over the 
range 275-295 nm (S 27 5_295) decreased weakly between the dry and 
wet period, with mean values of 0.0124 ± 0.0004 nm" 1 and 0.01 16 ± 
0.0003 nm" 1 , respectively. In contrast, S 3 5 0 _ 400 values tended to be 
higher during wet period (0.0152 ± 0.001 nm" 1 ) compared to dry 
period (0.0147 ± 0.0009 nm" 1 ). The difference in the S R ratio was 
however more marked between the two seasons, with lower values 
during the wet period (0.770 ± 0.049) compared to the dry period 
(0.845 ± 0.046). Note that the variations of DOM composition 
proxies were more marked when data are compared site to site. 
For example, the increase in a 350 in the Mbaere from 10.5 to 
35.0 m" 1 between dry and wet periods was related to clear 
decreases in the a 2 5o:a 36 5 ratio (from 4.3 to 3.7), S 27 5_295 (from 
0.0131 to 0.0117 nm" 1 ) and S R (from 0.860 to 0.799). Compared 
to the mainstem Oubangui River, 5 13 C values of DIC span a wide 
range of values and are typically more 13 C-depleted, with very low 
values in the Mbaere (—23.5 ± 2.0%o), intermediate values in the 
mainstem Lobaye ( — 16.1 ± 2.0%o) and in its minor tributaries 
( — 14.3 ± 2.3%o), and highest values in the Mbali and Mpoko 
River (—11.1 to — 10.5%o). 5 13 C DIC values are strongly negatively 
correlated with pC0 2 (Figure 7). With the exception of the Mpoko 
River, 5 13 C values of POC and DOC fall within the range expected for 
a dominance of C3 vegetation (-31.2 to -25.8%o for POC, -30.6 to 
—27.0 for DOC). The more turbid Mpoko River showed a 5 13 C PO c of 
— 22.8%o, indicating a substantial contribution from C4 vegetation in 
this subcatchment (no 5 13 C DO c data available for this site). 5 13 C 
values of POC and DOC pools were well correlated, and 5 13 C PO c 
values increased with TSM load (Figure 8). 

Calculated pC0 2 values for the different Oubangui tributaries 
were highly variable, ranging between 870 and 33500 ppm, being 
generally higher than those observed in the mainstem Oubangui 
River (Figure 9). Similarly, both dissolved CH 4 and N 2 0 concentra- 
tions were generally higher in tributaries than in the mainstem 
Oubangui River (Figure 9). 

Discussion 

Organic carbon. Rivers are often categorized in one of three types, 
mainly based on colour, distinguishing white-water rivers (alkaline 
to neutral pH, high sediment loads), black- water rivers (acidic, high 
DOC loads, very low sediment loads and dissolved ion 
concentrations), and clear- water rivers (low sediment loads but 
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Figure 3 | Seasonal variations of daily discharge (dotted lines) and (A) total alkalinity (TA) and the partial pressure of C0 2 (pC0 2 ), and (B) the carbon 
stable isotope signatures of dissolved inorganic carbon (5 13 C DIC ) in the Oubangui River at Bangui between March 2010 and March 2012. 



variable in terms of pH and alkalinity). This distinction is, however, 
not always straightforward (see e.g. 21). Consistent with the generally 
low mechanical erosion rates reported for the Congo Basin 11 , TSM 
concentrations in most of the tributaries were low. Only the Mpoko 
River which drains a wooded savannah region and some of the minor 
tributaries of the Lobaye, on the rainforest- savannah transition zone, 
showed slightly more elevated TSM concentrations (9.0-78.0 mg 
L" 1 ). While the Bodingue and Mbaere watersheds sampled here 
visually show characteristics of blackwater rivers, their DOC 
concentrations (1.5-7.1 mg L" 1 ) in the lower range of those found 
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Figure 4 | Seasonal variations of a 350 and the a 2 5o:a 3 65 ratio in the 
mainstem Oubangui River between March 2011 and April 2012. 



in the mainstem Oubangui (Figure 8), and generally lower than those 
reported for the mainstem Congo River (6.2 to 17.6 mg L _1 ; 13,19) or 
the Epulu River in the upper Congo Basin (5.2-9.0 mg L _1 ; Spencer 
et al. 2010), and lower than the range reported for blackwater rivers 
in the Amazon (7-40 mg L" 1 , 21). In all of the sampled tributaries, 
DOC remained the dominant form of organic C (OC), representing 
53-95% of the total OC pool. While sampling sites were not 
consistently similar between different sampling campaigns, DOC 
concentrations were generally higher during the 2012 wet season 
compared to the 2 dry season campaigns (e.g. for the mainstem 
Lobaye: 3.5 ± 0.9 mg L~\ n=5 during the wet season, 1.9 ± 

0. 1 mg L~\ n=5 during dry season; for the Mbaere: 6.4 ±1.0 mg 
L" 1 , n=2 during the wet season, 2.4 ± 0.8 mg L" 1 , n = 6 during dry 
season). Such higher DOC concentrations during wet conditions are 
in line with results from other sites in the Congo basin such as the 
mainstem Oubangui 1017 , mainstem Congo 1016 , and rivers in the Ituri 
Forest 14 . 

The ratio of spectral slopes S R and the 2L 25 o:a 365 ratio were 
developed as rapid and easy methods for characterizing cDOM. 
The S R ratio has been correlated to molecular weight and sources, 

1. e samples of greater allochthonous contribution with higher 
molecular weight DOM have lower S R values 22 . The a 2 5o:a 365 ratio 
has been correlated with the molecular size and DOM aromaticity 
with decreasing values relating to increasing molecular size 23 . In a 
recent study focusing on cDOM properties in 30 U.S Rivers, Spencer 
et al. 24 also showed that the a 2 5o:a 365 ratio was inversely related to the 
hydrophobic organic acid fraction (HPOA) of DOM which repre- 
sents the high molecular weight, aromatic- dominated fraction of 
DOM 25 . Therefore, low a 2 5o:a 36 5 ratio values (high HPOA fraction) 



SCIENTIFIC REPORTS | 4:5402 | DOI: 1 0.1 038/srep05402 



4 



o Oubangui River, 2010-2012 
• Oubangui tributaries 
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Figure 5 | Relationships between specific conductivity and (A) total 
alkalinity (TA), and (B) in situ pH for the mainstem Oubangui (grey 
circles; data from 2 years of monitoring), and tributaries of the Oubangui 
(black circles). 

indicate an allochthonous origin of DOM, whereas high a 2 5o:a 365 
values (low HPOA fraction) suggest an autochthonous sources (algal 
or microbial) or photodegraded DOM 26 " 27 . 

The ranges of 8275-295, S 350 _4oo and S R in the Oubangui system are 
comparable to that from other allochthonous dominated freshwater 
systems including temperate, Arctic, and tropical rivers 15 * 22 ' 28 - 30 . The 
range of the a 2 5o:a 36 5 ratio is comparable to that from other temperate 
and Arctic rivers 15 ' 28 . 

The large variations in stream S R and a 2 5o:a 36 5 ratios that occurred 
along with changes in DOC concentrations and water discharge 
during the high flow period in the Oubangui River and its tributaries 
clearly indicate that DOM pools are being mobilized that are differ- 
ent than those typical of base flow conditions (Supplementary Tables 
4 and 5). Similar changes in DOM optical properties across the 
hydrological cycle, especially in spectral slopes proxies, have also 
been reported in tropical 14 ' 29 and Arctic rivers 28 ' 30 , as well as in an 
agricultural watershed in California 31 . Highest DOC concentrations 
were linked with lower S R and a 2 5o:a 365 values, indicating an increase 
in the contribution of allochthonous sources to stream DOM due to 
greater surface runoff and leaching of organic rich layers during the 
high flow period 14 ' 30 . However, lowest DOC concentrations were 
linked with higher S R and a 2 5o:a 365 values, indicating that DOM is 
less aromatic in nature during the base flow period. This shift in 
DOM sources can be explained by the deepening of hydrologic flow 
paths through the soil profile and greater residences time of DOC in 
contact with soil subsurface microbial communities, leading to the 
mobilization of a more microbially processed DOM during base flow 
period 14 ' 30 . However, combining optical and fluorescence properties 
of cDOM, Yamashita et al. 29 proposed another explanation in the 
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Figure 6 | Relationship between dissolved organic carbon (DOC) 
concentrations and in situ pH for the mainstem Oubangui (grey circles; 
data from 2 years of monitoring), and tributaries of the Oubangui (black 
circles). 

interpretation of the increase of S R during base flow period reported 
in tropical rivers in Venezuela. These authors suggested that the 
increase of S R values could potentially result from a higher contri- 
bution of autochthonous source by enhanced plankton primary pro- 
ductivity during the low turbidity period. However, based on the 
previous study of Bouillon et al. 17 , such autochthonous source has 
not been identified as contributor of DOC in the Oubangui River, 
supporting the hypothesis of a deepening of hydrological flow path as 
driver of changes in DOM sources. 

The temporal variation observed in spectral slopes proxies in the 
Oubangui River are similar to those observed in the Epulu River 14 . 
During the wet period (April) the Epulu was characterized by lowest 
S R values (0.834 ± 0.028) close to those measured in the Oubangui 
River during high flow period occurring from September to 
December (0.832 ± 0.041). In the Epulu River, S R increased with 
decreasing discharge during post- and intermediary periods 
(November-February), with mean and maximal values of 0.962 ± 
0.056 and 1.066, respectively. Mean value of S R in the Oubangui River 
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Figure 7 | Relationships between the carbon stable isotope composition 
of DIC (5 13 C DIC ) and the partial pressure of C0 2 (pC0 2 ), for the 
mainstem Oubangui (grey circles; data from 2 years of monitoring), and 
tributaries of the Oubangui (black circles). 
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Figure 8 | Relationship between total suspended matter concentrations 
(TSM) and (a) 5 13 C signatures of particulate organic carbon (5 13 C PO c)> 
and (b)the contribution of particulate organic carbon (POC) to the TSM 
pool (%POC) for the mainstem Oubangui, tributaries of the Oubangui 
and samples collected throughout the Congo basin (literature data from 
Mariotti et al. 1991, Sigha-Nkamdjou et al. 1993, Coynel et al. 2005, 
Bouillon et al. 2012, Spencer et al. 2012). 



during base flow periods was higher, about 1.180 ± 0.116, with a 
maximal S R value of 1.323. S350-400 values showed the same temporal 
variations in both rivers, with lowest values in the Oubangui River 
compared to the Epulu River. The main difference between these two 
tropical rivers is the variations in S 27 5_295 values observed in the 
Oubangui River and to a lower extent in Oubangui tributaries. 
High flow period was associated with a clear decrease in 8275-295, 
whereas Spencer et al. 14 reported no seasonal variations in S 27 5_295 
across the hydrological cycle. 

Despite the similar temporal variations in DOC concentrations 
and DOM composition, spatial differences can be observed between 
the Oubangui River and its tributaries. The Oubangui River is char- 
acterized by a more important range of variations in whole cDOM 
parameters, including absorption coefficient, a 2 5o:a 365 (E2:E3), and 
spectral slope coefficients. This observation is quite different from 
those reported by Yamashita et al. 29 , where hydrological changes 
seem to affect cDOM properties in tributaries to a larger extent than 
the main river channels. 

S 13 Cpoc values showed a positive correlation with sediment loads 
(Figure 8), but were generally in the range expected for C3-domi- 
nated catchments. The 5 13 C values in the rivers draining the Ngotto 
Forest are relatively low ( — 32.2 to — 29.1%o), but vegetation in these 
systems can be expected to have low 5 13 C values given the effects of 
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Figure 9 | Boxplots of data on pC0 2 , CH 4 , and N 2 0 concentrations for 
the mainstem Oubangui River (data from 2 years of monitoring), and 
tributaries of the Oubangui (all sites and sampling seasons combined). 

pC0 2 data for the lower Congo River are from Wang et al. (2013). n 
indicates the number of data points. 

both high precipitation 32 and relatively closed- canopy structure 
(e.g., 33 ). The increasing pattern in 5 13 C PO c values with increased 
sediment loads (albeit still low on a global scale, 34) suggests higher 
sediment inputs in catchments where C4 vegetation is more substan- 
tial (i.e., from savanna landscapes). The negative relationship 
between TSM loads and %POC (the contribution of POC to the total 
particulate matter load, Figure 8) has been previously observed in a 
range of individual river basins (e.g., 35 ), and on a global scale 
(e.g., 4 ' 36 ). In our case, this correlation can be interpreted as reflecting 
the continuum from 2 contrasting end-members: (i) direct litter or 
organic-rich surface soil layers as observed in the majority of the 
rainforest rivers, and (ii) more soil- derived sediments as observed 
in the more turbid systems, with correspondingly lower %POC. 

Inorganic carbon. HC0 3 ~ in rivers is mainly derived from the 
weathering of carbonate and silicate rocks, with the relative 
proportions depending to a very large extent on the lithology of 
the drainage basin 37 " 38 . Gaillardet et al. (39, based on data from 8) 
estimated the relative contributions of carbonate and silicate 
weathering to the bicarbonate load of the Congo River and some 
of its tributaries (Oubangui, Sangha, and Kasai) and found a 
contribution of silicate weathering generally between 50-60% 
except for the Sangha where carbonate dissolution dominated 
(only 9% resulting from silicate weathering). Since we do not have 
data on Cl~ concentrations as required in their approach, we used the 
simple stoichiometric model of Garrels and Mackenzie 37 whereby the 
contribution to TA from carbonate weathering (TA carb ) is computed 
from Ca 2+ and Mg 2+ concentrations and the contribution to TA from 
silicate weathering (TA sil ) is computed independently from dissolved 
silicon (Si) concentrations, according to: 

TA carb = 2*([Ca 2+ ] + [Mg 2+ ] - [S0 4 2 "]) 
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Figure 10 | Relationship between (A) observed and modeled total 
alkalinity (TA) (see text for details), and (B) the carbon stable isotope 
composition of DIC (5 13 C DIC ) and the estimated contribution of TA 
derived from carbonate weathering to the observed TA. 

TA sil = [Si]/2 

S0 4 2 " allows to account for Ca 2+ originating from the dissolution 
of gypsum (CaS0 4 ), but this correction term was ignored in absence 
of S0 4 2 ~ data. In a companion study in the Tana basin (Kenya) we 
found S0 4 2 ~ to correspond to about 10% of Ca 2+ for a variety of 
streams and rivers 40 , which is likely an upper estimate for our study 
area given the more humid climate and the absence of evaporative 
environments. 

For the tributaries, there is a good linear regression between the 
modeled TA (TA carb + TA sil ) and observed TA (r 2 = 0.99, n = 36, 
Figure 10A). In the Oubangui mainstem, the modeled TA is well 
correlated to the observed TA but markedly deviates from the 1:1 
line as the values increase. The stoichiometric model of 37 assumes 
that the weathering of silicate rocks leads to a release of HC0 3 ~ and 
dissolved Si (DSi) according to a 1:2 ratio. This is true for common 
mineral forms such as olivine and albite, but the dissolution of some 
other mineral forms lead to a different HC0 3 ~:DSi ratio. For 
instance, the weathering of plagioclase feldspar (NaCaAl 3 Si 5 0 16 ) 
leads to the release of HC0 3 ~ and DSi in a 6:4 ratio. While no 
detailed information is available on the mineralogy of the underlying 
bedrock in our study catchments, they are underlain mostly by meta- 
morphic rocks in which feldspars are likely to be common consitu- 
ents. When we apply this ratio to the Oubangui mainstem data, the 
recomputed modeled TA in the high values range is in better agree- 



ment with observed TA, than based on the HC0 3 ~:DSi ratio of 1 : 2 
(Figure 10A). The contribution of carbonate rock weathering 
(%TA carb ) was estimated as the percentage of TA carb to total modeled 
TA (TA carb + TA si i), and ranged between 30% and 94%, encompass- 
ing the range of %TA carb computed by Gaillardet et al. 39 between 40% 
and 91%. %TA carb is positively correlated to TA indicating a lower 
contribution of carbonate rock weathering in the basins draining 
humid forest than savannah. 

Carbon in HC0 3 " originating from silicate rock weathering comes 
exclusively from C0 2 , while 1/2 of the C in HC0 3 ~ from carbonate 
rock comes from CaC0 3 and the other 1/2 from C0 2 . If the C0 2 
involved in the weathering comes from organic C degradation, 
S 13 C DIC should have a negative signature (as indicated by the 5 13 C 
signatures of DOC and POC), while marine CaC0 3 has a 5 13 C sig- 
nature close to 0%o 41 . 

This can explain the positive relationship between %TA carb and 
5 13 C DIC (Figure 10B). When applying a linear fit to these data 
(R 2 =0.69), the extrapolated 5 13 C DIC values for the end members 
where the contribution of carbonate weathering is 0 and 100% are 
— 31.3 and — 1.5%o, respectively, which corresponds well with the 
expected values for terrestrial vegetation in the rainforest biome 
and marine carbonates. 5 13 C DIC values as low as those found in some 
of the Ngotto forest rivers (as low as — 28.1%o) have, to the best of our 
knowledge, only been reported for a few rivers in the Amazon 
basin 42 " 43 , and can only occur under conditions when alkalinity is 
extremely low and hence where silicate weathering dominates. The 
strong correlation between 5 13 C DIC and pC0 2 (Figure 7) on the one 
hand, and the link between 5 13 C DIC and carbonate weathering dis- 
cussed above, suggest that the relative importance of silicate versus 
carbonate weathering may exert an important control on pC0 2 in the 
studied river systems. Such a mechanism is further suggested by the 
highly significant regressions (p<0.0001) between log(pC0 2 ) and 
the estimated contribution of carbonate weathering to TA 
(R 2 =0.38), or with Si/Ca 2+ ratios (R 2 =0.35), i.e.,rock weathering 
would explain close to 35% of the variance of pC0 2 . A recent analysis 
of factors controlling inorganic carbon speciation in North American 
rivers 44 concluded that while much of the variation in the spatial 
patterns of TA and pH can be explained by catchment processes 
and characteristics related to chemical weathering (e.g., precipita- 
tion, proportion of carbonate rocks), spatial variations in pC0 2 were 
mainly governed by in-river processes only indirectly related to the 
catchments (e.g. temperature through its effect on respiration). At 
the scale of our study, where factors such as temperature are much 
more homogeneous spatially and seasonally, it appears that the 
weathering regime could exert a partial control on the inorganic 
carbon speciation, and hence influence pC0 2 and water- atmosphere 
C0 2 exchange. It should be noted, however, that our calculated pC0 2 
data from rivers with low pH are likely to be overestimates, since a 
comparison of calculated pC0 2 data with those measured in the field 
with a direct headspace technique suggests that the former approach 
induces a substantial bias towards higher pC0 2 under conditions of 
low pH (<6.5) and/or very low TA (Abril et al., in preparation). 
While this compromises some of the absolute pC0 2 data from rivers 
in the Ngotto forest, it is unlikely to distort the overall patterns 
described above. 

Methane and nitrous oxide. The data compilation by Bastviken 
et al. 2 demonstrated the lack of CH 4 flux data from tropical 
systems, despite the fact that tropical freshwater systems are often 
claimed to show much higher CH 4 emissions than temperate or 
high-latitude systems. For the Oubangui tributaries, we have made 
no attempts to calculate diffusive CH 4 fluxes (due to the lack of data 
required to estimate gas exchange velocities, as outlined above for 
C0 2 ), yet it is evident that dissolved CH 4 concentrations are, on 
average, more than an order of magnitude higher than in the 
mainstem Oubangui River (Figure 9). Riverine CH 4 concentrations 
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in the few African catchments studied so far range drastically, with 
reported values between 1 and 6730 nM in the Athi-Galana-Sabaki 
River 45 , between 48 and 870 nM in the Comoe, Bia and Tanoe rivers 
in Ivory Coast 46 , and between 25 and 505 nM in the Tana River 
basin 47 . While CH 4 concentrations in the mainstem Oubangui 
show a pronounced seasonality (17, and Figure 2), a comparison of 
CH 4 data for the Lobaye and Mbaere between our dry and wet season 
campaigns shows no distinct differences. While the riverine network 
is clearly a consistent source of CH 4 to the atmosphere, the data 
suggest that diffusive CH 4 fluxes are highly variable in space and 
time and therefore difficult to upscale reliably. As reported earlier 
for the mainstem Oubangui River, dissolved N 2 0 levels were 
generally low (6.4-13.0 nM, except for the Kele, a small and 
shallow rainforest stream draining into the Mbaere where N 2 0 
concentrations reached 26.7 and 26.2 nM during the 2012 and 
2013 sampling campaigns, respectively). These values represent 
only a slight oversaturation with respect to atmospheric 
equilibrium (137 ± 22%, 380-402% for the Kele). A synthesis by 
Baulch et al. 48 found higher N 2 0 emissions to be associated with 
increasing N inputs, hence the low N 2 0 concentrations observed 
are in line with the pristine status of our study catchments and the 
observed low nitrate levels (median of 7.1 uM, range: 0-50 uM, n= 
32; own unpublished data). 

Concluding remarks. Overall, our results demonstrate that even 
within a relatively small set of subcatchments of the Oubangui 
River, biogeochemical characteristics of different rivers and 
streams show highly contrasting characteristics on the one hand, 
yet consistent co-variations between parameters on the other hand. 
According to Raymond et al. 3 and Regnier et al. 49 , one of the factors 
limiting our ability to better constrain the role of rivers in regional 
and global C budgets is that the current empirical database is too 
sparse to adequately resolve the diversity of soil types, inland waters, 
estuaries and coastal systems, and specifies the Congo Basin as one of 
the key areas of regional priority. The data collected and compiled 
here, while covering only a small subset of the Congo Basin, suggest 
that large data collection efforts will be required to adequately 
describe the aquatic biogeochemical variability for this large 
tropical basin. However, the systematic co-variation between some 
of the examined proxies of C origin and transformation, and the 
apparent link with certain catchment characteristics (lithology and 
vegetation) provide some important clues on the type of 
classification and characterization needed to upscale future data to 
the catchment scale. 

Methods 

The Oubangui River (Figure 1) is the second largest tributary of the Congo River 
(after the Kasai), with a length of 2400 km from the source (Uele River) to its 
confluence with the Congo River, and a drainage basin of 644000 km 2 , of which 
489000 km 2 (76%) is located upstream of Bangui 10 . The Oubangui catchment 
upstream of Bangui is dominated by dry wooded savannahs, with more humid forest 
situated downstream towards the confluence with the Congo mainstem. The mean 
annual precipitation in the catchment is —1400-1540 mm y" 1 10 . The hydrological 
cycle of the Oubangui is characterized by a single main flood peak and maximum 
discharge typically in October-November. Annual discharge has fluctuated between 
2120 m 3 s" 1 (1990) and 6110 m 3 s" 1 (1969). Although we did not find land-use data 
specifically for the Oubangui catchment, deforestation rates in the CAR are reported 
to be among the lowest in Africa (net deforestation of 0.06% y" 1 50 ). According to the 
Food and Agriculture Organisation (FAO) statistics (http://faostat.fao.org/), agri- 
cultural land use has remained stable between 1961 and 2010, and makes for —8% of 
the total land area. Considering a national population density of ~7 inhabitants 
km" 2 , these figures support classifying the Oubangui catchment as nearly pristine. 
Monitoring on the Oubangui in the city of Bangui (CAR, 4°21' N, 18° 34' E) was 
initiated in late March 2010, and was followed by approximately 2 -weekly sampling. 
A subset of data collected between March 2010 and March 201 1 were presented in 21 ; 
data on the subsequent year of monitoring (after which this programme ended) are 
included here. 

Tributaries of the Oubangui were sampled in March 2010 (dry season), March 
201 1 (dry season), and November 2012 (wet season). During the first field campaign, 
two sites were sampled along the Mbali River (Fig. 1), as well as the Mpoko River just 
south of Bangui (see also 13), and the main rivers draining the Ngotto Forest, Lobaye 



province (Fig. 1). The Ngotto Forest was revisited in 201 1 and 2012 for more extensive 
sampling on the main rivers (Mbaere, Lobaye, and Bodingue) and some of their 
minor tributaries. The Ngotto Forest, bordering on the Republic of Congo (a.k.a. 
Congo Brazzaville) and characterized by annual precipitation of — 1600 mm, is 
classified as a semi-deciduous rainforest and lies south of the wooded savannahs of 
the Guinea- Sudanian transition zone. Water levels were lower during the dry season 
surveys (March) than during the wet season (November), but water level differences 
are rather small (1.0-1.5 m), in line with the relatively stable discharge of the Lobaye 
River (monthly average flow ranging between 250 and 500 m 3 s-1 for the period 
1950-1984; data from the Global River Discharge Database (http://www.sage.wisc. 
edu/riverdata/, accessed February 2014). The Mbaere and Bodingue are the two main 
permanent rivers, which are bordered by extensive areas of seasonally flooded forests. 
The Bodingue River joins the Mbaere River prior to joining the Lobaye River, which 
lies on the northern side of the Ngotto Forest and represents the border towards drier, 
savannah landscapes to the North. The Lobaye River is an important tributary of the 
Oubangui River. The Ngotto Forest forms a protected area with high biodiversity, and 
very low human population densities 51 . The geology of the region is complex and 
includes igneous, metamorphic, and Mesoproterozoic- Neoproterozoic sedimentary 
formations that include calcareous sediments 52 . 

Water temperature, conductivity, dissolved oxygen (0 2 ) and pH were measured in 
situ with a YSI Pro Plus multimeter, whereby the 0 2 and pH probes were calibrated on 
each day of data collection using water saturated air and United States National 
Bureau of Standards buffer solutions (4 and 7), respectively. All sampling for larger 
rivers was performed from dugout canoes at —0.5 m below the water surface, in a few 
cases sampling was performed from a bridge (Mpoko) or from the river shore (Mbali). 
Samples for dissolved gases (CH 4 , N 2 0) and the stable isotope composition of dis- 
solved inorganic C (5 13 C DIC ) were collected with a Niskin bottle or a custom-made 
sampling bottle consisting of an inverted 1L polycarbonate bottle with the bottom 
removed, and —0.5 m of tubing attached in the screw cap 53 . 12 mL exetainer vials (for 
S 13 C DIC ) and 50 mL serum bottles (for CH 4 and N 2 0) were filled from water flowing 
from the outlet tubing, poisoned with HgCl2, and capped without headspace. 
Approximately 2000 mL of water were collected 0.5 m below the water surface for 
other particulate and dissolved variables, and filtration and sample preservation was 
performed in the field within 2 h of sampling. 

Samples for total suspended matter (TSM) were obtained by filtering 200-1000 mL 
of water on pre-combusted (4 h at 500°C) and pre-weighed glass fiber filters (47 mm 
GF/F, 0.7 um nominal pore size), and dried in ambient air during the fieldwork. 
Samples for determination of particulate organic C (POC), particulate nitrogen (PN) 
and C isotope composition of POC (§ 13 C PO c) were collected by filtering 50-300 mL 
of water on pre-combusted 25 mm GF/F filters (0.7 um nominal pore size). The 
filtrate from the TSM filtrations was further filtered on 0.2 um polyethersulfone 
syringe filters (Sartorius, 16532-Q) for total alkalinity (TA), DOC and o 13 C DOC (8- 
40 mL glass vials with Polytetrafluoroethylene coated septa). TA was analysed by 
automated electro-titration on 50 mL samples with 0.1 mol L" 1 HC1 as titrant 
(reproducibility estimated as typically better than ± 3 umol kg" 1 based on replicate 
analyses). The partial pressure of C0 2 (pC0 2 ) and DIC concentrations were com- 
puted from pH and TA measurements using thermodynamic constants of Millero 54 as 
implemented in the C02SYS software 55 . It should be noted, however, that this 
approach is expected to overestimate pC0 2 in some of the blackwater rivers (see 
Discussion). 

For the monitoring data on the Oubangui River, C0 2 exchange with the atmo- 
sphere was calculated as F= kaApC02 where k is the gas transfer velocity, oc is the 
solubility coefficient for C02, and ApC02 represents the difference in partial pres- 
sure of C0 2 between water and air. Monthly averages corresponding to time of 
sampling of atmospheric pC0 2 from Mount Kenya (Kenya, -0.05 °N 37.80°E) were 
retrieved GLOBALVIEW-C02 database (Carbon Cycle Greenhouse Gases Group of 
the National Oceanic and Atmospheric Administration, Earth System Research 
Laboratory), k values were calculated using two approaches, one based on the rela- 
tionship between wind speed and k600 (i.e., k normalized to a constant temperature of 
20°C), and an alternative approach using a gas transfer parameterization based on 
depth and water current that was developed for rivers and streams; details can be 
found in Bouillon et al. 17 . No C0 2 flux estimates were made for the tributary sampling 
sites, since both approaches would be hampered by a lack of adequate ancillary 
information - data on river depth and flow velocity are lacking, and the varying degree 
of canopy closure would compromise the use of modeled wind speed data. 

For the analysis of 5 13 C DIC , a 2 ml helium (He) headspace was created, and H 3 P0 4 
was added to convert all DIC species to C0 2 . After overnight equilibration, part of the 
headspace was injected into the He stream of an elemental analyser - isotope ratio 
mass spectrometer (EA-IRMS, ThermoFinnigan Flash HT and ThermoFinnigan 
DeltaV Advantage) for 5 13 C measurements. The obtained 5 13 C data were corrected 
for the isotopic equilibration between gaseous and dissolved C0 2 as described in 
Gillikin and Bouillon 56 and measurements were calibrated with certified reference 
materials LSVEC and either NBS-19 or IAEA-CO- 1. Reproducibility of 5 13 C mea- 
surements was 0.2%o or better. Concentrations of CH 4 and N 2 0 were determined via 
the headspace equilibration technique (20 mL N 2 headspace in 50 mL serum bottles) 
and measured by gas chromatography (GC, 63) with flame ionization detection (GC- 
FID) and electron capture detection (GC-ECD) with a SRI 8610C GC-FID-ECD 
calibrated with CH 4 :C0 2 :N 2 0:N 2 mixtures (Air Liquide Belgium) of 1, 10 and 
30 ppm CH 4 and of 0.2, 2.0 and 6.0 ppm N 2 0, and using the solubility coefficients of 
Yamamoto et al. 57 for CH 4 and Weiss and Price 58 for N 2 0. 

25 mm filters for POC, PN and 5 13 C PO c were decarbonated with HC1 fumes for 
4 h, re-dried and packed in Ag cups. POC, PN, and § 13 C P oc were determined on the 
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abovementioned EA-IRMS using the thermal conductivity detector (TCD) signal of 
the EA to quantify POC and PN, and by monitoring m/z 44, 45, and 46 on the IRMS. 
An internally calibrated acetanilide and sucrose (IAEA-C6) were used to calibrate the 
5 13 Cpoc data and quantify POC and PN, after taking filter blanks into account. 
Reproducibility of 5 13 C PO c measurements was better than ±0.2%o. Samples for DOC 
and 5 13 C DO c were analysed either on a Thermo HiperTOC-IRMS, or with an 
Auroral030 TOC analyser (OI Analytical) coupled to a Delta V Advantage IRMS. 
Typical reproducibility observed in duplicate samples was in most cases <±5% for 
DOC, and ±0.2%o for 5 13 C DO c- For a subset of samples from the Ngotto Forest, we 
compared DOC concentrations and § 13 C data for samples filtered on 0.2 um (as 
described above) and only on 0.7 urn GF/F filters; for both variables no significant 
differences were found (n= 15, two-tailed paired t-test, in both cases p<0.0001 at 95% 
confidence interval, data in Supplementary Table 1). 

Samples for major element concentrations and Si were filtered through a 0.45 um 
polyethersulfone (PES) or polycarbonate (PC) filter, and measured by ICP-AES (Iris 
Advantage, Thermo) or ICP-MS (Perkin Elmer Elan 6100). 

Samples for the analyses of cDOM spectral characteristics were prepared similarly 
as for DOC analyses, but without H 3 P0 4 addition. cDOM samples were stored in 
amber glass vials with PTFE-coated septa. Absorbance measurements were recorded 
on a Perkin- Elmer UV/Vis 650S using a 1 cm quartz cuvette. Absorbance spectra 
were measured from 190 to 900 nm at 1 nm increments and noise instrument was 
assessed measuring Milli-Q water as blank. The correction for scattering, index of 
refraction and blank was performed by fitting the absorption spectra to the data over 
the range 200-700 nm according to the following equation: 

A^Aoe-^-^+K, (1) 

where A^ and A 0 are the absorbance measured at defined wavelength X (nm) and at 
reference wavelength X 0 =375 nm, respectively, S is the spectral slope (nm -1 ) that 
describes the approximate exponential decline in absorption with increasing wave- 
length, and K is a background offset. According to Johannessen and Miller 59 , the offset 
value was then subtracted from the whole spectrum. The fit was not used for any 
purpose other than to provide an offset value. After correction, absorption coeffi- 
cients were calculated according to the relation: 

a x = 2.303 xAx/L, (2) 

where aX is the absorption coefficient (m _1 ) at wavelength X, A is the absorbance 
corrected at wavelength X, and L the path length of the optical cell in meters (0.01 m). 

Several optical indices were calculated to investigate CDOM properties, including 
spectral slopes over the ranges 275-295 nm (8275-295) and 350-400 nm (S35o_4oo)> the 
slope ratio (SR) and the a 2 5o:a 36 5 ratio. Spectral slopes were determined using linear 
regression for the log-transformed a spectra for the intervals 275-295 nm and 350- 
400 nm and the slope ratio SR was calculated as the ratio of S 27 5_295 to S 350 _4oo 22 . The 
a 2 5o:a 3 65 ratio (also called E2:E3 ratio) was calculating using absorption coefficients at 
the appropriate wavelengths. 

Literature data on specific parameters were obtained either from tabulated values 
in the relevant papers, data supplied by authors, or digitized from figures using 
PlotDigitizer v2.6.1 (http://plotdigitizer.sourceforge.net/). 
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